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Abstract 

Arguments in favor of a large magnitude ( at least two- three 
times bigger than its factorized value) of the mixed vacuum condensate 
(qG®G®q) are given. The analysis is based on the strict inequalities 
which follow from the QCD sum rules method and on very plausi- 
ble phenomenological assumptions like rag s > for the few lowest 
exited Qq states in a heavy quark limit thq — > oo. The same argu- 
ments show the suppression of the SU(3) symmetry breaking effects 
for vacuum condensates when the additional gluon fields are included. 



1 On leave of absence from Budker Institute of Nuclear Physics, 
Novosibirsk, 630090, Russia. 

e-mail addresses:arz@mail. physics. smu.edu, ariel@sscvxl.ssc.gov 



1. Formulation of the problem 

It is widely believed that the QCD vacuum has a very complex structure. 
To understand this structure is a very ambitious goal which assumes the solu- 
tion of a whole spectrum of tightly connected problems, such as confinement, 
chiral symmetry breaking phenomenon, the U(l) problem, the 9 dependence 
problem and many, many others. 

A less ambitious purpose is the phenomenological study of the QCD 
vacuum structure from the known experimental data. The most appropri- 
ate analytical approach which makes contact between fundamental theory, 
QCD, from the one side and observable variables from the other side, is the 
QCD sum rules [Q,0]. This method has been used extensively for about 
fifteen years for various purposes, mainly for the extraction of information 
on hadronic properties. It is important to stress that the enormous wealth 
of data referring to the low energy hadronic physics obtained in this way, is 
determined by only a few fundamental characteristics 

(qq) ~ -(0.25GW) 3 
(^G 2 ^) ^ 1.2 ■ lO^GeV (1) 

7T 

A a 

(qiga^G^—q) ~ m 2 (qq), m 2 ~ 0.8GeV 2 

Here q describes any of the light quark fields u, d, s and G^ u is the field of 
colored gluons. 

The value of the chiral condensate (qq) has been known for a long time 
H and the value of the gluon condensate has been first extracted from the 
analysis on charmonium system The value of the mixed quark-gluon 
condensate has been extracted first from the analysis of the nucleon [|J. 

It is clear that the QCD vacuum is a very complicated system and a 
few global dimensional parameters (H) give only some insight on the scale 
of vacuum fluctuations, but do not at all exhaust all characteristics of the 
vacuum structure. 

For studying of the fine aspects of the theory it is necessary to know 
the averages of more complicated operators like ((qq) 2 ), (qG^G^q), (G 3 ), 
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(G*), etc. Comparing these averages with global characteristics ([!]) gives a 
quantitative information of the granulated structure of the vacuum, charac- 
terizes the correlation between vacuum fluctuations of large (~ R con f) and 
small (~ GeV~ x ) size, gives some insight on the role of the interaction of the 
various vacuum fluctuations, etc. For instance, if the vacuum would be ho- 
mogeneous one, we would get more or less an uniform distribution for gluonic 
G^y and quark qq vacuum fields. The factorization hypothesis should work in 
this case perfectly well and predicts that (G* v ) ~ (G%,) , ((qq) 2 ) — (qq) 2 '■ At 
the same time, some granulated distribution of the vacuum fluctuations, like 
instantons, gives a very strong deviation from this prediction. In particular, 
the instanton liquid model predicts the violation of the factorization 

hypothesis on the level of the factor 3-10 (numerical factor depends on color 
and Lorentz structure of the operator). 

Let us note that the status of this hypothesis for four-quark condensates 
((qTq) 2 ) (dimension six) has been examined in many different ways (see || 
for review). Several channels have been studied with the special task of 
detecting deviations from the factorization for the four-quark operators. The 
conclusion is as follows: in vector and axial- vector cases ( T — j^, 7^75 ) 
deviation from factorization formula cannot significantly exceed ~ 20%. On 
the contrary, it can be shown |7] that this situation is not universal and in 
exotic cases with (r = <y^ u , 75, 1) , a noticeable violation of factorization does 
occuif| 

The purpose of this letter is to obtain phenomenological, model indepen- 
dent information on mixed vacuum condensates of dimension seven (qG^ u G^ v q) , 
which are minimal possible operators where factorization procedure can be 
applied ((qG^ v G^ v q) ~ (qq) ■ (G^G^)) and checked. In addition, we study 
the SU(3) breaking effects for these mixed condensates. Some motivation 
(apart from the academic interest concerning QCD vacuum structure) for 
the considering these VEVs will be explained in the conclusion. 

Before we proceed, let us describe briefly the main idea of the approach 
which we follow in this paper. As mentioned above, the QCD sum rules make 
contact between vacuum characteristics through their condensate values and 
some hadronic properties. It turns out that the heavy-light (Qq) quark sys- 
tem is a very sensitive to the high dimensional condensates like (qG ^q) , 
(qG^G^q) ... Thus, some hadronic characteristics are determined by these 

2 Somc qualitative arguments on why it should be so, are given in the last section. 
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vacuum expectation values (VEVs). 

Once this is said, the standard line of reasoning (very roughly) is the 
following. The magnitude of the condensate gives, on one hand, the char- 
acteristic scale where deviation from the asymptotic behavior is still under 
control. On the other hand the same scale defines some dimensional hadronic 
parameters which, in principle, can be extracted by applying some fitting 
procedure. 

We go in the opposite direction. We are not interested in calculation of 
any hadronic parameters in this paper. We are not calculating, in particular, 
rriB s or rriB u in the limit itlq oo. We are not going to prove the power of 
the method by demonstrating that mB s —mB u — lOOMeV in agreement with 
experimental data. Instead, we assume that the QCD sum rules describe 
this system well enough. In particular they should explain the inequality 
m Bs > m Bu which we assume is still valid in a heavy quark limit mq — > oo. 
Because of the sensitivity of the corresponding sum rules (which we are going 
to consider) to the higher dimensional condensates, this assumption (mg s > 
mB„) gives a very nontrivial restriction on the absolute values of these VEVs. 
Such restrictions are the main result of this paper. 

The next, part of the paper is devoted to the description of the method 
we are going to use. We consider the well-known p — <fi system and find 
the inequality which follows from the fact that > m p . The obtained 
restriction is a very weak and not interesting. The only justification for 
including this section to the paper is to formulate an idea with the simplest 
example. The third part of the paper is its main part where we apply the 
previously formulated idea to Qq system. Some comments are given in the 
conclusion. 

2. The Method. 

We introduce the jff* = S'j^s current with the meson quantum numbers 
and define the corresponding polarization operator in a standard way: 

i J e**dx{T{J<f>(x), JW(0)}> = {q,q u - q 2 g^P W (Q 2 ), Q 2 = -q 2 .(2) 

The calculation of the correlator for Q 2 — > oo with the power corrections 
taken into account is done in the standard way, and after the Borelization 



3 



procedure the sum rule takes the following form M: 

1 r s T . , M 2 f1 6m 2 87r 2 (m s ss) 

- / e m* ImP w ( s )ds = — r{l f H x f 7 

W V ; 4tt 21 M 2 M 4 



+ 3 M 4 81 s M6 W 

The matrix elements entering to eq. (|3|) are discussed above and we assumed 
the factorization for four-quark operator with ( Y = 7 M , 7^75 ) [0. In 
principle, we could analyze this sum rule in a standard way in order to 
reproduce the well-known experimental data for meson. It turns out that 
it is very useful to consider a slightly different sum rule by taking the first 
derivative with respect to M 2 . In this case, the left hand side of eq.(BI) is 

S 

proportional to - / e ~TZ*ImPW>{s)sds. If the meson would saturate these 
sum rules, we would get: 

2 - J e~M^ ImP^(s)sds ^{1+power corrections} 
m <t> — \ „ s_ ; , w ; . = m2 ri , - 7- (4) 



Je~t%ImPW(s)ds ^{l+power corrections} 



The same procedure can be done for the p meson with substitution m s — > 
m U} d and (ss) — > (mm) |. Now it is time to introduce the new parameter 
i?(M 2 ) defined as 

, 9N i fe~M^I m P^(s)sds ± [ e'th ImP^(s)ds 
R(M ) = 7tJ _ g — g _ g — — (5) 

i/e M^ImP^>(s)ds ±fe M s ImP(p\s)sds 

If the lowest state would saturate the corresponding sum rule, we would get 

2 

i?(M 2 ) ~ ^ > 1. (6) 

The important thing for us in this paper is that i? > 1. Intuitively it is 
clear and related to the fact that the mass of the strange meson (or meson 
built from the strange quarks) is bigger than the mass of its nonstrange 



3 we assume an exact isotopical invariance and do not distinguish (uu) and (dd). 
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partner. More important, this inequality R > 1 is still correct even without 
the assumption on saturating the dispersion integral by the lowest state. 

We have to make a few very mild assumptions about the first 2 — 3 lowest 
states, whose contribution to the dispersion integral is not suppressed by the 
weight factor e - *^. For those states we assume that the masses and the 
corresponding coupling constants of the strange particles bigger than their 
non-strange analogies ( for example, f K i > f^, > f p i, m^i > m pl , i = 
1, 2, 3...). An enormous wealth of experimental data and physics intuition tell 
us that it is definitely should be true. In this case, we are quite sure that the 
inequality R > 1 holds even when saturation by the lowest resonance does 
not take place. 

Taking into account the previous discussing and substituting the theo- 
retical parts of the corresponding sum rules into the expression fl5|) we get 

& 

R = l + AR 1 + AR 2 (7) 



m 2 AT ^ 9 (m s ss) 448 , (uu) 2 — (ss) 2 . . 

We divided AR = AR 1 + AR 2 on purpose. Ai? x in the above equation de- 
notes the perturbative contribution to AR and AR 2 denotes the correspond- 
ing nonperturbative part. They have absolutely different physical meaning 
and have different origin. It is clear that AR\ is related to the perturbative 
contribution. The insertion (in order to calculate the polarization operator 
(HD of the nonzero quark mass (m s ) to the asymptotic quark loop clearly 
decreases the phase volume. It automatically leads to the increasing of the 
corresponding hadronic mass. The fact ARi > is the explicit demonstra- 
tion of it. The contributions to ARi come exclusively from the small distance 
physics. In a sense, AR\ describes some trivial kinematical factors; there is 
no deep, nonperturbative physics, involved with ARi. 

In the rest of this paper, we will be interested in the nonperturbative 
part of the decomposition (0), namely AR 2 . One may wonder whether the 
decomposition (|7|) is unique. We refer to the original paper |L| for a detail 

4 We assume that the SU(3) breaking effects are small enough in order to use an ex- 
pansion with respect to this small parameter. It is clear that in the exact SU(3) limit, the 
ratio R = 1. 
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discussion of Wilson operator expansion (OPE) in given context, but here 
we want to stress that the decomposition (|7]) has the same status as the 
separation of all field fluctuations in large and small scales within an (OPE) 
and so it is well defined and based on solid ground. 

There is no general theorem concerning the sign of Ai? 2 |]. However, the 
physical picture of the origin of the hadron mass assumes its nonperturbative 
nature. Let us note, in passing, that in the chiral limit, the mass of hadrons 
comes exclusively from nonperturbative physics. 

Thus, it is natural to assume that the same, nonperturbative effects are 
responsible for the mass difference of a hadron and its partner with strange 
quarks. It is parametrically true in the chiral limit m s — > for the p — 
system (|3|). In this case, A_Ri ~ m 2 s , AR 2 ~ m s ^> ARi. As we know m\ — 
m 2 p ~ m s and so, the mass difference in this case is determined exclusively 
by the nonperturbative part. 

In Qq system we are going to discuss in the next section, both parts 
Ai?i,Ai?2 have the same order of magnitude ~ m s . However, we believe 
that the nonperturbative effects, as well as perturbative ones, work in the 
same direction in order to split the strange and nonstrange partners. All 
sum rules, I am aware of, satisfy this assumption, and the formal expression 
which is given by formula 

AR 2 > 0. (9) 
For — p system (|8|) this restriction leads to the following inequality: 

\ m s (ss) \> ^^[(uu) 2 - (ss)% (10) 

where Mq ~ ~ lGeV 2 is the characteristic scale when the corresponding 
sum rule does make sense. Let me briefly explain this point. As usual, the 
analysis of the sum rules is somewhat of an art. We could take parame- 
ter Mo to be a very large; in this case the next power corrections could be 
definitely neglected, but the resulting inequality would be trivial and non- 
informative. If we take M to be too small, there is no quarantee that a 
higher power corrections will not dominate. Our choice of M , when the 
main power correction is less then 20% of the asymptotically leading term, is 

5 The sign of ARi is always positive by kinematics. 
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the safe one in order to neglect other power corrections; usually this parame- 
ter numerically is very close to the position of the lowest resonance in a given 
channel. In particular, for the p meson, Mq (p) ~ 0.6GeV 2 , for if* and <j) it is 
M 2 (K*) ~ 0.8GW 2 , M o 2 (0) ~ l.OGeV 2 respectively. Starting exactly from 
these points, the "theoretical" and "phenomenological" parts of the sum rules 
fit each other. From the other side, we are still in the transition region, and 
very sensitive to VEVs. To be safe enough for the numerical estimations, 
we will choose the parameter Mq to be bigger than any lowest resonance in- 
volved in our analysis. In particular, for p — <p system we take Mq ~ lGeV 2 . 
For this choice we are quite sure that the next power corrections are small 
enough and under control. 

Now, a few remarks regarding ex.(|It]). First of all, as it was expected, 
both parts of this expression are proportional to m s . In addition, because 
we know the SU(3) breaking effects for the condensate (see reviews 
[§] and references on previous papers therein), 

(mm) — (ss) 

A = _ , V ' ~ 0.2, 11 

(mm) 

and one can see that the formula ([LT]) is trivially satisfied. No new results can 
be obtained in this particular case. However, in the next section we will derive 
an analogous inequality where instead of (mm) — (ss), there enters the higher 
dimensional operators like (uiga^G^u) — (siga^G^s) or (ug 2 G fMU G^ u u) — 
(sg 2 G lJiU G^ l/ s) which are less familiar. In this case we can study the influence 
of the gluonic fields on the SU (3) breaking parameters and check the standard 
conjecture about weak influence of the gluonic fields on the SU(3) breaking 
parameters, i.e. 

A=?=A 9 =?=A 99 =?=... (12) 

A = (uiga^G^u) - (siga^G^s) _ (ug 2 G^ u G^ u u) - (sg 2 G ilv G lJLV s) 

9 (uiga^Gfivu) ' 99 (^G^G^u) 

If we knew the A g and A gg parameters independently, the expression analo- 
gous to ( P~D|) would give us information on absolute value of the condensate. 
In the example under consideration it leads to 

m M 2 27 

(0.25GeV) 3 ~\ (mm) |< . r „ ~ (0.55GeV) 3 . (13) 
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We expect an accuracy of such inequalities is about ~ 20%, as the standard 
accuracy of the sum rules approach. The next section is devoted to Qq system 
where some more interesting operators (instead of (uu)) will enter into the 
game. In this case an inequality analogous to (|1^) will be less trivial and 
much more interesting. 

3. Heavy-light quark system 

Introduce the J = Qijsq current with the quantum numbers of the pseu- 
doscalar meson and define the corresponding polarization operator in a stan- 
dard way: 

i J e iqx dx(T{J(x) + , J(0)}) = P{q 2 ). (14) 
We use below the technique proposed in pi, i.e. the energy e is used instead 



of q 2 : q 2 = (Mq + e) 2 , e <C Mq. The advantage of this technique is its 
simplicity; the disadvantage is the existence of large 1/Mq corrections for the 
real D, 5-mesons. However, as we mentioned, we are not going to extract 
any information about real hadrons in this paper, we are going to study 
vacuum of QCD. In this case, with assumption that no qualitative changes 
will occur in the limit tuq — > oo, it is much easier to live in this imaginary 
world. Besides that, let us note, that the first calculations of fjj, fs within 
the sum rule approach have been done using exactly this technique |10 1 , [|i~l 



independently the same system has been investigated by the standard method 
in ref. [0]. In addition, this technique is very useful to make contact with 



heavy-quark effective theory |13fl , see recent paper [JT4 . 

A few remarks regarding new variable e instead of q 2 : as usual, we are 
going to calculate the correlation function in the nonphysical region, e < (it 
is analogous to q 2 < in the standard method). In the limit E = — e — > oo, 
(but E <C Mq in order to use heavy quark approximation), we can use the 
standard OPE and Borel transformation, defined by operator B. In given 
case B takes the following form: 

B= lim , 1 ,, £"(-- ^-)" (15) 

E,n->oo-A=M=const.(n-l)\ dE' 

where the energy E is positive and measured from the threshold. The physi- 
cal meaning of the Borel transformation in this case much more clear than in 
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the relativistic case. Indeed, when we are passing from E to M we actually 
go from Green function at given energy to the time Green function. The con- 
tribution of the each state with energy E{ to the dispersion integral is given 

_£i 1 

by formula e m and -g plays the role of time it. It is clear that the OPE 
expansion ™ at M — > oo corresponds to t — > 0, i.e. asymptotically small 
times. To study a low lying resonance we have to calculate the correlator at 
small enough M, which corresponds to large distance physics (large times t). 
It is easy to check the following properties of operator B: 

b(^) = jkh^iw^ B{EHnE) = H)* 1 *^ ( 16 ) 



a 1 1 

B = — exp m 

E + Ei M 

After these few general remarks let us calculate the nonperturbative, 
Borel transformed part of P np (M) with strange s quark in the current J = 

Qi"f 5 s : 

1 i i 13 - 2!M 3 

Pnv = \ — (ss) H (m s ss) H n(siq<J uu G uu s)}, P v = ' (1 7) 

p 2M AM 4M 2 M M /J p 2M tt 2 

We removed the common factor ^jj on purpose- the contribution of any 
physical resonance to the dispersion integral has the same factor and 
after Borelization takes the following form: 

ftl f ImP(s)ds f2( Ml 1 1 ^ 

B n J s + Q* fA ^ ] 2~MM- Q e ' (18) 

where we used the standard definition for f P 

M 2 

(0\Qt7 5 q\P Q ) = f PJ f. (19) 
Q 

Besides that, the factor i in (18) is related to Borel transformation (|15|), and 
comes from nonrelativistic 5 function :5(s — M P ) = 5(2Mq(E — E R )) = 

j^5(E — Er). We included to eq . (fP7D the formula for perturbative, asymp- 
totically leading part as well. The reason to do so just to make sure that 
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at the scale we are going to discuss, the power corrections are much smaller 
than the leading term. But essentially, the only information we need to know 
about this term is its dimension, which can be found by pure dimensional 
arguments. The inequalities we are interested in and which follow from the 
our main nonperturbative assumption @ do not depend on the constant in 
front of P p and we will not discuss the perturbative contribution any more. 

Let us note that P(M) does not depend on Mq. From this observation it 
is easy to reproduce the well known result fpMq ~ 1 of the effective heavy 
quark theory. 

Now we are ready to repeat all steps of the previous section which lead 
us to the eq. (|10D from the definition @and assumption (||). For the given 
correlator (JTJ|) the analogous formula looks as follows 

1 5 

(ss - uu) - 7^( m sS s ) ~ -^-p^iga^G^s - uiga^G^u) > 0, (20) 
which we would like to rewrite in the following form: 

5mn fn s x 

-A < — + MA, (21) 

12M 9 3 ' y ! 

where Wq ~ Q.RGeV 2 is parameter for mixed condensate ([!]); m s ~ Q.l'bGeV 



and A ~ 0.2 flTT|). Before we proceed with numerical estimations, let me 
make a few comments on (^). It is clear, that as M — > oo this inequality 
is trivially satisfied because (ss — uu) = A\(qq)\ > 0. It becomes less trivial 
in the transition region. In this case the last term, proportional to A 9 , 
contributes with the opposite sign Q The requirement of positivity of AR2 
impose some limitation on the absolute value of A g . So, our main task is to 
find parameter Mq where higher power corrections can be safely neglected, 
but inequality will be still interesting. As we said above, we are not going 
to extend the sum rules approach on the region below the first resonance. 
In given case, from the QCD sum rules analysis in the heavy quark limit 
(see recent papers Jl4]] , |L7|] , [16|]) and from experimental data (we expect 
that B system already very close to heavy quark limit), we can learn that 
M Bu — Mq ~ 0.5GeV and M Bs — Mq ~ O.QGeV. So, our choice 

M ~ 0.6GeV (22) 



3 There is no doubt in sign A g p5| 
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should be considered as safe enough. Power corrections, as can be checked 
from (|i"7f) already small. With this remark in mind we will get the following 
(rather weak, again) restriction: 

A g < 0.3. (23) 

It is definitely in agreement with what we knew before H15| , |T8[| . Moreover, 
from the first estimate |H| , we could suspect that the SU (3) breaking effects 



for mixed condensate are much weaker than for the standard (qq) condensate. 
It would be very interesting to check whether it is indeed true by some 
independent analysis. In order to do so, we want to consider the correlators, 
which include derivatives. The reason for this choice will be clear soon. 

Introduce the J M = Qij^iD^q current with the quantum numbers of the 
pseudoscalar meson and define the corresponding polarization operator in a 
standard way: 

i J e^dx(T{J^x) + , J(0)}> = q»P {1 \q 2 ). (24) 

Here = <9 M — igA^^- is the covariant derivative acting on the light quark. 
The physical sense of the corresponding matrix element 

M 2 

(OlQi^iD^PQ^)) = q J P -f( Xq ). (25) 

M Q 

can be understood from the observation that (x q ) determines the mean mo- 
mentum fraction carried by the light q quark in the meson Pq with total 
momentum q^. This information is very important for description of the 
asymptotic behavior of exclusive processes; we refer to original paper |TlJ on 



this subject, but now we want to remark that (x q ) ~ <C 1 and so, the 
corresponding common factor will occur in phenomenological as well as 
in the theoretical parts. It leads, of course, to the increasing dimension of 
the P^(q 2 ) on one unit in comparison with P(q 2 ) from (|i~4l). 

As was mentioned in |TT] the main power correction to the theoretical 
part of (p4|) is determined by the mixed operator in the chiral limit. So the 
corresponding correlation function is very sensitive to this operator and we 
have a good chance to improve our limit (p3|). This is the main motivation 
to consider (^4]) in this paper. 
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By repeating again all steps which lead us to the formula fllTD , we get the 
following Borel transformed expression for PW; 

7-.m 1 r 4a s M _£o 1 _ 1 „ . 

p i = 2mmq { ~^t ( ' s)(1 - e ' M) + 4 {ms ' s) + m { * l9a » vG » vS) ~ 



m s ,.. n n „i 1 3-3!M 4 

32^ (S8 ^ G - S)} ' P W = 2MMq • (26) 

Few comments are in order. First of all, as was expected, the leading at 
M — > oo term proportional to chiral condensate (qq) is suppressed by the loop 
correction ~ — and numerically small. This fact increases sensitivity to the 
next VEV (siga^G^s) . As the second remark, let us note the appearance 
of the continuum threshold parameter so in the nonperturbative part P^j} . 
This is the standard parameter which accompany all sum rules, but usually it 
occurs only in the perturbative parts. Its physical meaning can be seen from 
the standard model for the spectral density as the resonance plus continuum 
which starting at sq. The reason of appearing such parameter in our case is 
the high dimensions of our currents. The consequence is the growth of the 
spectral density ImP$9(s) with s. From the QCD sum rules analysis in the 



heavy quark limit |T^|, p7|1 , Ul6[ , we can learn that 



So = m R + 0.5GeV. (27) 

It is interesting to note that this "mnemonic rule" works amusingly well in 
all known massive channels. Anyhow, because the condensate contribution 
comes with the small factor ~ — , the influence of so is negligible. As usual 
we included to the formula ( [ZED the perturbative part as well in order to 
demonstrate the smallness of the power correction at the scale Mq — O.GGeV. 

Let us rewrite the final inequality, analogous to (|20|) and which follows 
from our main conjecture @: 

(1 — e m — e m ){ss — uu) — (m s ss) — 

vr 3M n 



5 3/77- 

■isiga^G^s - uiga^G^u) + lA sigo^G^s) > 0. (28) 
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It can be rewritten in much more clear way, using our notations for A and 



8M 9 s 7T V 3M ; 16M 2 ' V ; 

For numerical estimations we use the standard set of parameters already nor- 
malized at the low normalization point: /i ~ lGeV; m s ~ 0.15Ge^; m 2 , ~ 
0.8GeV 2 ; a s (/x) ~ 0.34. Exactly these parameters have been used in the 



analysis of qQ system in the heavy quark limit |L6| with the special attention 
on the loop corrections and we prefer not to change it. With these remarks 
in mind we will get the following very strong restriction on A g : 

A g < 0.15. (30) 

This is in agreement with our feeling gathered from the absolutely inde- 



pendent consideration |15j| . The relation (|30|) means that the gluon inser- 
tion into the chiral condensate makes (siga^G^s) numerically very close to 
(uiga^G^u). 

Few words on the accuracy. We expect that the next power corrections 
can not exceed 10 — 20% anyhow. Besides that, we keep only linear SU(3) 
breaking terms ~ m s . In particular, is treated as A and not as 

A(l + A). The term ~ A 2 is order m 2 and should be discarded (or we must 
gather all these corrections at once). We expect that these next corrections 
~ m 2 can not exceed 10 — 20%. Finally, the reason for our belief in ( |30| ) is 
that the analysis [|n| of the quite different system of Goldstone bosons, gives 
the same conclusion. 

As a last remark concerning formula (^), we want to stress one more 
time that this is not a sum rule in the usual sense. We are not fitting it 
in order to extract some information on hadron properties. Instead, we do 
believe that the corresponding sum rules reproduce mass spectrum correctly, 
in particular, strange hadron heavier than its non-strange partner. It leads 
to some relations between VEVs. We expect that the obtained inequalities 
should hold in the region of variable M where sum rules are suppose to be 
working. 

With this remark in mind and in order to get some feeling on the high 
dimension condensate (qg 2 G fJiU G^ u q) , let us introduce the current J^ u = 
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Qi^iD ^iD ' u q and consider polarization operator 

i J e iqx dx(T{J^(x) + , J(0)}> = g M g,P (2) (g 2 ) + g^P±. (31) 

We will be interested in consideration of q^q v kinematical structure only. The 
corresponding matrix element 

_ M 2 
^\Qi^D^D v q\P Q {q ll ))^q v q ll jp-^{x 2 )+g lJ , v ... . (32) 

Mq 

determines the mean square of the longitudinal momentum carried by light 
quark. It is clear that (x 2 ) ~ Mr and so we expect an additional suppression 

q M Q 

of the correlator in comparison with the previous case and increasing 
of its dimension on one unit more. An explicit calculation supports this 
expectation and the corresponding expressions for the perturbative and 
nonperturbative P^ parts of (^TJ) take the form 

r-.m 1 r 20L a M . _Hl 47T — 2 

^(m.Hg^^) - ^p(0,)h = 2 J M , 4 ' . (33) 

In this formula the operator (O q ) has dimension seven and defined as follows 

(O q ) = {qg 2 G^ u G^ v q) - ^{qg 2 a^G^a^, u ,G^ u i) (34) 

where as usual G^ v = G"^. The factorized value of this condensate is given 
by 

1 1 

6 + Yl 



(O q ) f = (g 2 G%G%) ■ (qq) • (« + ^) (35) 



where the coefficient | comes from the factorization of the first term (^) and 
Y2 comes from the last one. Let us note that the sign of (O) is determined 
uniquely from the inequality obtained below and independently from fl9| ; it 



coincides with its factorized value (|35|). Thus we define 

(O) = (0) f ■ K, K>0 (36) 
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and our task now is to find some restrictions on the coefficient of nonfactor- 
izability K. To do so we follow the standard procedure which bring us to 
eqns.d^^): 

7?r 2 ,« s ^2 v „ a . 32Ma 8 , 1 



BO • K • A 99 > -m s ■ ^(1 - e~* - £Ur*)+ 



Q M 2 > ^ MW 99 5?r V 3M 

m s m§ 64vra s A 

i»r--w l<w>| - (37) 

With the standard input parameters it gives: 

K-A gg > 0.38. (38) 

As we argued above, the insertion of an additional gluon field into the con- 
densate (qq), leads to the suppression of the SU(3) symmetry breaking ef- 
fects in the obtained mixed condensate (q~Gq) of approximately one quarter: 
A g /A ~ 0.75. We expect that the same suppression holds for each of the 
next gluon insertions. Thus, for 0.75A g < A gg < A g we obtain: 

K > 2.5(forA ga = A g ); K > 3A(forA 9g = 0.75A g ), (39) 

which is our main result. The independent analysis of the different system 
(Goldstone bosons, not Qq mesons) supports this result and demonstrates 
the violation of factorization by a factor 3. 

4. Conclusion 

Thus we have found two qualitative phenomena from the analysis of the 
Qq system (very sensitive to vacuum structure) in the heavy quark limit 
M Q -> oo: 

• SU(3) breaking effects are being suppressed with the insertion of gluon 
fields into condensate. In particular, we expect that the following link of 
inequalities is correct A > A g > A gg > ... I do not know whether any QCD 
vacuum models can (at least) qualitatively describe this behavior. 

• The factorization hypothesis does not work for mixed, dimension seven 
operator (qGGq). It is not a big surprise, because we faced the analogous 
phenomenon early for the four-fermion condensates with an "exotic"[] 
Lorentz structure. 



7 This term has been borrowed from the Shifman Introduction to the book g] 
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In conclusion, I want to give some qualitative remark which elucidates 
the reason why some VEVs must violate factorization. First of all, let me 
remind you that in some special cases the QCD sum rules do not work in 
principle. In particular, in the pseudoscalar quark channel 

z iqx dx(T{J^\x), J W (0)}) = P (lT \Q 2 ), Q 2 = -q 2 , J = uij 5 d (40) 

the standard sum rule can not reproduce the residue (0|J|7r) = f n ■ 2GeV 
which is known exactly from PCAC, and which has much bigger scale than 
QCD sum rule can provide. The reason for that became clear very soon after 
QCD sum rules were invented and related to the fact that the limited dy- 
namical information encoded by means of the different VEVs is not sufficient 
for describing the vacuum ± channels. So called direct instanton contribu- 
tions ]2T1 play a decisive role in these channels. Let us call these currents as 
"exotic" ones, to emphasis their difference in comparison with the "normal" 
currents. The scale provided by QCD sum rules analysis in the "normal" 
cases does describe all residues correctly. 

Now let us consider any nondiagonal correlator, constructed from the 
"normal" and "exotic" currents. On the one hand, there are no direct 
instanton contributions; on the other hand, the "exotic" matrix element 
(0| J\tt) = f^2GeV is several times larger than the normal scale f n -m p , which 
the sum rule can provide. Thus, the information about big scale connected 
with contributions of direct instantons, clearly filters through into ImP(s) 
via the residue r n ~ (\J exotlc \7r) {7i\J normal \) . The question arises, in which 
way can one guarantee a large scale of the quantity rv, if the corresponding 
sum rules do not allow direct instanton contribution? From our point of view 
the answer is that some VEVs which enter into the corresponding correlator 
become numerically large. Precisely in such a way ( through enhancement of 
the vacuum condensates) the information of large value of the scale filters 
through into correlators, which do not allow direct instanton contributions. 

• Few words on applications obtained results. First of all, as was ex- 
plained in the text, these mixed VEVs determine the momentum distribution 
between q and Q quarks in the heavy-light system. In any calculation with 
qQ system involved, this information is very useful, see e.g. 21] . 



As a second remark, I would like to note that the transverse quark distri- 
bution in the light mesons is determined mainly by high-dimensional conden- 



sates studied above, see |0| for a more details. As is known, the transverse 
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size dependence plays a key role in color transperency physics ||22|| . 
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